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ATMOSPHERIC  CORRECTION  OF  AVIRIS  DATA  IN  OCEAN  WATERS 

Gregory  Tcrric  and  Robert  Arnonc 
Naval  Research  Laboratory,  Stennis  Space  Center 

1.  INTRODUCTION 

Hypcrspectral  data  offers  unique  capabilities  for  characterizing  the  ocean 
environment.  The  spectral  characterization  of  the  composition  of  ocean  waters  can 
be  organized  into  biological  and  terrigenous  components.  Biological  photosynthctic 
pigments  in  ocean  waters  have  unique  spectral  ocean  color  signatures  which  can  be 
associated  with  different  biological  species  (Arnone  et  al,  1986).  Additionally, 
suspended  sediment  has  different  scattering  coefficients  which  result  in  ocean  color 
signatures.  Measuring  the  spatial  distributions  of  these  components  in  the  maritime 
environments  provides  important  tools  for  understanding  and  monitoring  the  ocean 
environment.  These  tools  have  significant  applications  in  pollution,  carbon  cycle, 
current  and  water  mass  detection,  location  of  fronts  and  eddies,  sewage  discharge 
and  fate  etc. 

Ocean  color  has  been  used  from  satellite  for  describing  the  spatial 
variability  of  chlorophyll,  water  clarity  (k^^),  suspended  sediment  concentration, 
currents  etc.  Additionally,  with  improved  atmospheric  correction  methods,  ocean 
color  results  have  produced  global  products  of  spectral  water  leaving  radiance  (L,,) 
(Feldman  1989).  Ocean  color  results  have  clearly  indicated  strong  applications  for 
characterizing  the  spatial  and  ten.poral  variability  of  bio-optical  oceanography 
(Arnone  et  al,  1991). 

These  studies  have  been  largely  the  results  of  advanced  atmospheric 
correction  techniques  applied  to  multispectral  imagery.  The  atmosphere  contributes 
approximately  80%  -  90%  of  the  satellite  received  radiance  in  the  blue-green 
portion  of  the  spectrum.  In  deep  ocean  waters,  mtodmum  transmission  of  visible 
radiance  is  achieved  at  490nm.  Conversely,  nearly  all  of  the  light  is  absorbed  by  the 
water  at  wavelengths  greater  than  about  650nm  and  thus  appears  black.  These 
spectral  ocean  properties  are  exploited  by  algorithms  developed  for  the  at.-nospheric 
correction  used  in  satellite  ocean  color  processing. 

The  objective  of  this  study  was  to  apply  atmospheric  correction  techniques 
that  were  used  for  processing  satellite  Coastal  21one  Color  Scanner  (CZCS)  data  to 
AVIRIS  data.  Quantitative  measures  of  from  AVIRIS  will  be  compared  wth 
ship  ground  truth  data  and  input  into  bio-optical  models. 

2.  METHODS 

Atmospheric  correction  of  CZCS  assumes  a  linear  summation  of  the  total 
radiance,  L,(>i),  sensed  at  the  sensor  (Gordon  et  at,  1983,  Gordon  and  Clark,  1979). 
The  atmospheric  contribution  is  dividea  into  the  Rayieigh,  L,(A),  component 
(molecular  scattering)  and  the  aerosol,  L,(A),  component  (Mie  scattering).  The 
volume  scattering  of  the  water,  L„(A)  (ocean  color)  is  transmitted  through  the  air 
sea  interface,  t(A).  All  of  these  parameters  arc  related  to  the  total  radiance,  L,(A) 
by 

L,(X)  =Z,,a)  (1) 
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Assuming  that  the  signal  at  670nm  is  due  entirely  to  atmospheric  effects  and  the 
aerosol  optical  thickness  r,(670)  is  proportional  to  r,(X),  equation  (1)  can  be 
wriuen  in  terms  of  L^(A),  at  nadir,  by 

r£JA)=L,(X)-i,(A)--^^(^)V,(670)-L/670)).  (2) 


'.'.here  Fj  is  the  solar  flux  and  n  is  the  angstrom  coefficient. 

This  relationship  was  applied  to  AViRlS  data  by  expanding  this 
uoiieralized  nadir  equation  to  account  for  the  angular  atmospheric  path  radiance 
irom  the  sensor  altitude.  The  rayleigh  component  L,(A)  was  computed  using  a 
multiple  scattering  model  (Gordon  et  al,  1988).  Both  solar  and  sensor  zenith  and 
azimuth  angles  were  computed  for  each  pixel  in  the  scene  according  to  the  date, 
time  and  location  of  data  collection.  The  solar  fluxes  were  assumed  constants  and 
the  angstrom  coefficient  was  varied. 

Four  spectral  channels  were  used  in  processing  CZCS  ocean  color  data 
which  were  centered  at  443nm,  520nm,  550nm,  and  670nm  with  a  band\»ndlh  of  20 
nm.  The  443nm  channel  represents  the  chlorophyll  absorption  wavelength  which  is 
inversely  proportional  to  concentration  (Clark  1980).  The  520nm  and  550nm 
channels  characterize  the  hinge  point  in  bio-optical  ocean  waters  and  arc  related  to 
sediment  scattering.  These  bands  were  simulated  using  AVlRlS  bands  by  averaging 
the  two  adjacent  bands  for  each  wavelength.  Additionally,  a  750nm  channel, 
analogous  to  the  CZCS  750nm  band,  was  used  as  the  atmospheric  reference  channel 
because  of  a  problem  with  bottom  reflectance  at  670nm,  especially  in  the  shallow 
coastal  regions.  Each  of  the  channels  were  calibrated  for  spectral  radiance  and  for 
vignetting  effects  using  the  supplied  correction  data. 

The  atmospheric  model  was  applied  to  an  AVIRIS  scene  collected  off 
Tampa  Bay  on  April  3,  1990.  This  image  has  many  different  ocean  regimes 
characterizing,  shallow  coastal  waters  with  a  uniform  bottom  reflectance,  suspended 
sediment  in  the  water,  absorbing  "yellow  substance",  and  variability  of  coastal 
phytoplankton  and  clear  waters.  Ship  measurements  were  collected  by  the 
University  of  South  Florida  (Carder  et  al,  1991)  coincident  with  the  overflight  at 
three  stations  in  which  L,(A),  chlorophyll,  and  the  diffuse  attenuation  coefficient 
were  measured. 

3.  RESULTS 

The  443nm  channel  is  severely  contaminated  by  the  atmosphere  and  sensor 
noise  such  that  little  ocean  features  arc  observed.  Atmospheric  scattering 
diminishes  in  the  520nm  and  550nm  channels  and  ocean  features  can  be  recognized 
although  not  quantitative.  Atmospheric  correction  using  the  670nm  channel  as  the 
reference  was  applied  to  the  scene  ranging  from  clear  waters  to  the  coastal 
environments.  In  the  near  shore  coastal  waters,  the  AVIRIS  data  was  shown  to  be 
correlated  with  the  bathymetry  since  the  reflectivity  of  the  white  bottom  clearly 
dominated  the  water  leaving  radiance  signal.  The  610nm  channel  was  shown  to  be 
correlated  with  bottom  reflectivity  and  therefore  is  not  a  good  channel  for  only 
characterizing  the  atmosphere,  Bccau.se  of  these  non  zc'c  ''’TOrsiu  components, 
these  is  an  overcorrcction  for  the  atmosphere.  This  results  in  negative  radiance  in 
areas  where  the  bottom  reflectance  dontinates  the  670nm  channel.  Since  this  is  not 
possible  the  images  were  processed  and  negative  values  clipped  to  zero.  The 
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Assuming  that  the  signal  at  670nm  is  due  entirely  to  atmospheric  effects  and  the 
aerosol  optical  thickness  »’,(670)  is  proportional  to  r,(A),  equation  (1)  can  be 
written  in  terms  of  Lv(A),  at  nadir,  by 

rL„(X)=L,(X)-L(A)--^^(^)''{L,(670)-L,(670)).  (2) 

where  Fj  is  the  solar  flux  and  n  is  the  angstrom  coefficient. 

This  relationship  was  applied  to  AVIRIS  data  by  expanding  this 
generalized  nadir  equation  to  account  for  the  angular  atmospheric  path  radiance 
irom  the  sensor  altitude.  The  rayleigh  component  L,(A)  was  computed  using  a 
multiple  scattering  model  (Gordon  et  al,  1988).  Both  solar  and  sensor  zenith  and 
azimuth  angles  were  computed  for  each  pixel  in  the  scene  according  to  the  date, 
time  and  location  of  data  collection.  The  solar  fluxes  were  assumed  constants  and 
the  angstrom  coefficient  was  varied.  i 

Four  spectral  channels  were  used  in  processing  CZCS  ocean  color  data 
which  were  centered  at  443nm,  520nm,  550nm,  and  670nm  with  a  bandwidth  of  20 
am.  The  443nm  channel  represents  the  chlorophyll  absorption  wavelength  which  is 
inversely  proportional  to  concentration  (Clark  1980).  The  520nm  and  550nm 
channels  characterize  the  hinge  point  in  bio-optical  ocean  waters  and  arc  related  to 
sediment  scattering.  These  bands  were  simulated  using  AVIRIS  bands  by  averaging 
the  two  adjacent  bands  for  each  wavelength.  Additionally,  a  750nm  channel, 
analogous  to  the  CZCS  750nm  band,  was  used  as  the  atmospheric  reference  channel 
because  of  a  problem  with  bottom  reflectance  at  670nm,  especially  in  the  shallow 
coastal  regions.  Each  of  the  channels  were  calibrated  for  spectral  radiance  and  for 
vignetting  effects  using  the  supplied  correction  data. 

The  atmospheric  model  was  applied  to  an  AVIRIS  scene  collected  off 
Tampa  Bay  on  April  3,  1990.  This  image  has  many  different  ocean  regimes 
characterizing,  shallow  coastal  waters  with  a  uniform  bottom  reflectance,  suspended 
sediment  in  the  water,  absorbing  "yellow  substance",  and  variability  of  coastal 
phytoplankton  and  clear  waters.  Ship  measurements  were  collected  by  the 
University  of  South  Florida  (Carder  et  al,  1991)  coincident  with  the  overflight  at 
three  stations  in  which  L„(A),  chlorophyll,  and  the  diffuse  attenuation  coefficient 
were  measured. 

3.  RESULTS 

The  443nm  channel  is  severely  contaminated  by  the  atmosphere  and  sensor 
noise  such  that  little  ocean  features  are  ob.scrvcd.  Atmospheric  scattering 
diminishes  in  the  520nm  and  550nm  channels  and  ocean  features  can  be  recognized 
although  not  quantitative.  Atmospheric  correction  using  the  670nm  channel  as  the 
reference  was  applied  to  the  scene  ranging  from  clear  waters  to  the  coa.stal 
environments.  In  the  near  shore  coastal  waters,  the  AVIRIS  data  was  shown  to  be 
correlated  with  the  bathymetry  since  the  reflectivity  of  the  white  bottom  clearly 
dominated  the  water  leaving  radiance  signal.  The  670nm  channel  was  shown  to  be 
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chaiacterizing  the  atmosphere.  Because  of  these  non-zero  G70nm  components, 
these  is  an  overcorrection  for  the  atmosphere.  This  results  in  negative  radiance  in 
areas  where  'he  bottom  reflectance  dominates  the  670nm  channel.  Since  this  is  not 
possible  the  images  were  processed  and  negative  values  clipped  to  zero.  The 


